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ABSTRACT 



A one-dimensional mixed-layer model of ocean thermal 
structure (Camp, 1976) was first modified to determine the 
effects of a salinity profile on density structure and 
secondly to observe the added effects of surface salinity 
flux. The model was tested in August and November-December 
1974 with (1) hydrocast data at Ocean Stations HOTEL and 
PAPA and (2) salinity profiles statistically derived from 
historical data for six ocean stations and adjusted to 
correspond to an observed BT. Inclusion of salinity struc- 
ture did not significantly affect mixed-layer depth or 
temperature predictions during the summer period. During 
the winter period, salinity structure tended to inhibit 
deepening yielding a slightly warmer, shallower mixed 
layer. Surface salinity flux altered significantly the 
thickness of the isothermal layer, with decreases in thick- 
ness during precipitation (E-P<0) and a tendency for in- 
creased thickness during some periods of net downward heat 
flux when E-P>0 . 
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I. 



INTRODUCTION 



It has been established that the upper layer of the 
ocean responds to atmospheric forcing on the synoptic and 
diurnal time scales [Elsberry and Camp (1978), Garwood 
(1977)]. The success of those models which incorporate 
only mechanical mixing, convective overturning, and vertical 
diffusion in predicting the time evolution of the mixed 
layer (Camp and Elsberry, 1978) has shown the short term 
ocean boundary layer response to be predominately one-dimen- 
sional in many ocean regimes . Although the general model 
[see Camp (1976)] contains both temperature and salinity 
effects, Camp and Elsberry (1978) considered density as a 
function of temperature only, and assumed buoyancy flux to 
the atmosphere to be attributable to the heat flux alone. 
Johnson (1977) followed with an evaluation of the model's 
ability to forecast ocean thermal structure based on atmos- 
pheric forcing derived from synoptic-scale data. Again, 
due to lack of salinity data, Johnson neglected salinity 
effects and therefore had to adjust initial temperature pro- 
files to remove positive gradients. It is the purpose of 
this research to study salinity effects on the upper ocean 
thermal structure using the Camp (1976) model by making com- 
parisons with the results of Johnson (1977). 

Miller (1976) included salinity in the mixed layer model 
of Kraus and Turner (1967). Using BOMEX Period III upper 
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ocean soundings, he showed, for periods up to 10 days, the 
effects of salinity on mixed-layer temperatures and depths 
and the effects of precipitation on surface salinity and 
temperature in the tropical Atlantic. He concluded that 
the most significant effect of inclusion of salinity was 
the reduced cooling rate of the mixed layer due to reduced 
entrainment of water from the pycnocline into the base of 
the layer. This is a result of a larger density jump at 
the base of the layer when a salinity jump corresponds to 
a temperature jump than would be the case if density were 
a function of temperature alone. This type of density 
structure can also lead to a positive temperature gradient 
(temperature increasing with depth) immediately below the 
mixed layer (i.e. the mixed-layer temperature is cooler 
than that of the underlying water). Miller also concluded 
that the effects of precipitation were greatest under condi 
tions of light winds and heavy rainfall during which a 
shallow stable layer is produced at the ocean surface. 

Considering Miller's findings, it is apparent that 
salinity is important in determining the density structure 
when applying a mixed-layer model to the entire ocean. Fur 
thermore , although the relative importance of salinity 
effects on the short term evolution of the density profile 
may not be significant, one would anticipate that a long 
term imbalance between evaporation and precipitation contri 
butes to the buoyancy flux and could lead to significant 
density changes. Additionally, sound speed in the ocean 
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and salinity variations 



varies by 1.2 to 1.3 ms ^ (°/oo) ^ 
in the near-surface layers can thus have significant ef- 
fects on the sound-speed profile. These effects would be 
greatest in regions of varying surface salinity flux and 
advection (e.g., doldrums, subpolar regions). 

The importance of salinity in the model and the paucity 
of concurrent temperature and salinity data, even at ocean 
weather stations, raise the question of how to incorporate 
salinity effects into the initial conditions and to then 
verify the results when the only available information is 
a temperature profile from a bathythermograph. It was the 
objective of this study to test the salinity contribution 
to the density structure observed based on two different 
approaches: (1) utilization of available hydrocasts when 

both temperature and salinity profiles are available, and 
(2) utilization of bathythermographs and typical tempera- 
ture/salinity profiles derived from the historical data 
files and adjusted to be compatible with an observed BT. 
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II. MODEL MODIFICATIONS 



Camp (1976) and Miller (1976) presented the governing 
equations for the vertical redistribution of heat and turbu- 
lent kinetic energy in the upper ocean. These equations, as 
applied to a well-mixed layer, were simplified by appropri- 
ate assumptions and closure was achieved by parameterization 
of the dissipation and shear production terms of the turbu- 
lent kinetic energy budget. At the surface, the boundary 
conditions were specified in terms of the wind stress, inso- 
lation, sensible and latent heat fluxes, back radiation, and 
evaporation minus precipitation. The boundary conditions at 
the base of the layer are then computed by the model in 
terms of the fluxes of momentum, heat, and salt. In develop- 
ing the numerical scheme from the expressions, Camp neglected 
salinity effects because of a lack of salinity data. For 
this study, the model was modified to account for the contri- 
butions of salinity to the density structure and the 
buoyancy flux. 

A linear equation of state was assumed, which, when ex- 
pressed in terms of buoyancy, B, is 

B = agT - BgS (2-1) 



where 



a 



1_ 9_p 

p 9T 
o 



and 
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1_ _9_p 

p 9S 
o 
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are taken as constant; a : 2.5 x 10 (C) and 3 = 

7.5 x 10 -4 (°/oo ) _1 after Camp (1976). The assumption that 
3 is constant is quite good. However, density is not a 
linear function of temperature and a can be taken as con- 
stant only over a limited range of temperatures. The den- 
sity profile is represented in the model as temperature and 
salinity pairs to a depth of -NAZ by N isothermal and iso- 
haline slabs of AZ(= 2.5 m). 

The potential energy (PE) per unit is then expressed as 



PE = 




(aT 



3S) ZdZ 



( 2-2 ) 



where D is a fixed depth greater than the maximum penetra- 
tion of the vertical turbulent processes . Assuming advec- 
tive processes and the variation due to compressibility are 
negligible, the change in potential energy resulting from 
mixing the layer from a depth -NAZ to a depth -(N+1)AZ is 



APE(N) = h P o gN(AZ) 2 [a(T N -T N+1 )-3(S N -S M+1 )] . (2-3) 

If a(T^-T^ + ^ )-B ( ) > 0, the water column is stable 
and APE (N) is the amount of turbulent kinetic energy re- 
quired to mix the layer to a depth -(N+1)AZ. If a(T^-T^ T+ ,.) 

- B(S^-Sj^ + ^) < 0, the column is unstable and will overturn, 
generating turbulent kinetic energy by free convection (re- 
presented as APE (N)). 

c 
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In each one hour time step the model first accounts 
for surface fluxes, and the resulting temperature profile 
and surface (N=l) salinity are determined by 



T (t+At) = T 1 (t)+ -^[-Q T (0,t*) + Q s (0,t*)-Q (-AZ,t*)] 

o p 

T (t+At) = T (t)+-^[Q s ( -NAZ,t* ) - Q s ( -(N+ 1) AZ,t*)] 

o p 

S 1 (t+At) = S 1 (t)+ S 1 (t) ||[E(t*) - P(t*)] 



(2-4) 



where Q is solar radiation, Q T is the sum of latent, sen- 

O x 

sible, and back radiation, E is the evaporation rate, and 
P is the precipitation rate during the interval t* = At. 
Insolation is distributed vertically in the water column 
by absorbing 50% in the first meter and absorbing the remain- 
ing 50% as exp(-yZ). The extinction coefficient, Y, was 
assumed constant at 0.002 cm - '*" after Johnson ( 1977 ). 

If the resulting buoyancy profile is unstable, the water 
column is mixed according to 



,(i + l) 
N+l 



( NT ^ i) + T^)/(N+1) 



(i + 1) 
N + l 



(Ns[ l} + S^^ 1) /(N + 1) 



(2-5) 



until the profile is stable. The temperature and salinity 
profiles are then isothermal and isohaline to a depth 
-(N+1)AZ. The potential energy generated by this free con- 
vection is found from 



APE 

c 




-D 



NN-1 

(a6T- B6S ) ZdZs» ^ 

N= 1 



APE (N) 
c 



( 2 - 6 ) 
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where 5T and SS are the changes in temperature and salinity 
resulting from mixing and the depth of free convection is 
(NN)(AZ) < D. 

If the resulting buoyancy profile is stable, additional 
kinetic energy must be expended to further mix the water 
column. The total amount of energy, E^,, available for mix- 
ing the first N levels with the N+l level is 



E™ = E + E - E 
T m c p 



(2-7) 



where E^ is the mechanically generated turbulent kinetic 
energy, E c is the convectively generated turbulent kinetic 
energy, and E^ is the amount of energy previously expended 
to mix the water column to depth -NAZ. 

E is calculated by 
m J 



E (N) = [p w* exp ( -NAZ/H) ] At 
m o- r 



( 2 - 8 ) 



with w... = average friction velocity in the ocean over the 
time interval At and H = scale depth of 50 m. This depth- 
dependent formulation follows the approach of Elsberry, 
Fraim, and Trapnell (1976). Setting a scale depth of H=°° 
would correspond to the original Kraus and Turner (1967) 
model in which the net generation minus dissipation was in- 
dependent of the depth of the layer. 

E is considered to be the fraction of APE available 
c c 

for entrainment and is calculated as 



NN-1 

E (N) = -R V APE (N) (2-9) 

c / j c 

N=1 
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where R = 0.15 following Camp (1976) and Johnson (1977), 
and represents the percentage of convectively generated 



energy not dissipated by viscous forces, 
from 



N-l 

E (N) = APE (i) 

p Z— / m 

i = NN 



E is calculated 
P 



( 2 - 10 ) 



If E,p > APE^ for N levels, then there is energy avail- 
able to mix to level N+l. If E_ < APE , the level is par- 

I m 

tially mixed following Thompson (1976) and Camp (1976). 
Thus , let 

ne t (n) 

a = (N+1)APE (N) 
m 



Q = [aQ N+1 - (N-a) Q N ]/N (2-11) 

Q'= a Q N + d-a) Q n+1 



and set = Q for i=l,2,...,N and Q^ + -^ = Q' where Q is tem- 
perature and salinity. Note that the mixing algorithms re- 
quire that the isothermal and isohaline layers be the same 
depth. The mixed-layer depth is determined from the tempera- 
ture profile only by defining it as the depth at which the 
temperature is 0.2°C less than T^ , where T^ has been assumed 
to be equal to T , the sea-surface temperature. Light wind 
conditions, strong heating, or heavy precipitation could re- 
sult in some difference between T^ and T . The pycnocline 
is assumed to be weakly diffusive and the change of tempera- 
ture and salinity below the well-mixed layer is specified by 
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8T 

at 





( 2-12 ) 



_as 

at 



where the vertical eddy diffusion coefficient, A , is taken 



value is consistent with many observations. It is much 
greater than the values of molecular conductivity and dif- 
fusion. The processes involved are not well understood but 
may be due to breaking of internal waves, double diffusion, 
and/or lateral advection on a small scale. 



2 - 1 

to be 0.5 cm -sec after Haney and Davis (1976). This 
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III. DATA 



This study used atmospheric and bathythermograph data 
extracted from the Fleet Numerical Weather Central (FNWC) 
historical data file by Johnson (1977). Additionally, the 
12-hour accumulated precipitation field produced by FNWC's 
primitive equation model was extracted. Oceanographic sta- 
tion data for Ocean Station HOTEL were obtained from the 
National Ocean Data Center (NODC) and for Ocean Station PAPA 
from the Institute of Ocean Sciences (1975). The data at 
OS HOTEL consisted of 18 Nansen cast observations and the 
data at OS PAPA consisted of three hydrographic observations 
and 11 STP ( salinity- temperature-pressure ) observations for 
the time frames indicated in Table I. Additionally, daily 
surface salinity samples were available at OS PAPA for the 
period 06 December to 31 December 1974. 

To take advantage of the more numerous bathythermograph 
(BT) observations, it .was proposed to use some type of 
historical salinity profile in conjunction with a BT to 
initialize the model. One possibility was to use a profile 
generated from climatological means at different depths. 

This proved unfeasible in that the profiles contained very 
little detailed information on the salinity structure in 
the upper levels of the ocean. Another possibility, using 
standard temperature-salinity relationships for a water mass, 
was rejected since these relationships fail in the boundary 
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layer because of changes caused by the surface fluxes. The 
third possibility was to utilize the Hydroclimatological 
Data Base Retrieval (HYDAT) system installed at FNWC (Ocean 
Data Systems, Inc., 1977). HYDAT retrieves hydroclimatologi- 
cal data from the historical and synoptic data bases. FNWC ' s 
historical files contain expendable bathythermograph and 
Nansen cast data. The synoptic file contains the most re- 
cent 72 hours of BT reports. HYDAT then statistically ana- 
lyses these data as described in Appendix A and selects a 
temperature profile that is considered typical for the time 
period (>_ 1 month) specified. 

If the typical profile is a Nansen cast, the salinity 
profile is output with the temperature profile. If the 
typical profile is a bathythermograph, then no salinity pro- 
file is given. In this case, it may still be possible to 
obtain a salinity profile by searching manually through the 
HYDAT output for a Nansen cast sounding whose sea-surface 
temperature, deep temperature, and heat content come closest 
to fitting the median values. For this study, the parameters 
were weighted such that the fit of sea-surface temperature 
and heat content were considered more important than the 
fit of deep temperature. In some locations, Nansen casts 
were sparse and their parameters deviated from the median 
values. Thus, a minimum requirement for considering a Nan- 
sen cast as a typical profile was that the three parameters 
must fall within the range between the first and third 
quartiles. Proceeding in this manner, salinity profiles 
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were obtained for all six ocean weather stations in December 
and for three ocean weather stations in August (Table I). 
These salinity profiles were assumed to be typical for these 
ocean station locations and time of year. More ideally, one 
would like to have a data bank of hydrocast data and an ex- 
tracting system which would take a bathythermograph and 
search the historical data for the closest fitting tempera- 
ture profile and then use the corresponding historical 
salinity profile as the typical profile. 

Once the typical profiles were obtained, the question 
arose as to how to utilize them. The initial observed tem- 
perature profile likely has structural details that should 
be reflected in the initial salinity profile. These details 
include the depth of the well-mixed layer and positive tem- 
perature gradients. Thus the typical salinity profile should 
be mixed as deep as or deeper than the BT and the salinity 
gradient should be. such that positive temperature gradients 
do not lead to an unstable density profile. The typical pro- 
files for a month were therefore treated in two ways to 
generate a salinity profile that would correspond to any 
specified BT for that month. 

The first method was to combine the typical temperature 
and salinity profiles, T^(Z) and S^CZ) into a typical buoy- 
ancy profile B t (Z), by 

B t (Z) = agT T (Z) - 6gS T (Z) (3-1) 
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and the buoyancy profile was checked for stability. Then 
the mixed-layer depth, the depth where the temperature is 
0.2°C less than the surface temperature, of the typical tem- 
perature profile and of the bathythermograph were determined. 
If the layer depth of the T,p(Z) profile was less than that 
of the bathythermograph, the buoyancy profile was mixed 
according to equations (2-5) to the layer depth of the 3T. 
This mixing is necessary because of the assumption in the 
model that the isothermal and isohaline layers are mixed to 
the same depth. If the mixing were not done, the model 
would use a certain amount of input energy in mixing the 
salinity profile to the level of the temperature profile, 
whereas, in reality, it is likely that the actual salinity 
profile was already mixed to the deeper level. If the layer 
depth of the T^,(Z) profile was greater than that of the BT, 
no mixing was done. If a positive temperature gradient was 
encountered in the BT before the layer depth was reached, 
the BT was considered to have a mixed-layer depth equal to 
zero. Finally, with any given bathythermograph observation, 
T q (Z), and the assumed buoyancy profile B,p(Z), a correspond- 
ing salinity profile, S , was found from 

agT (Z) - B t (Z) 

S (Z) = ^ ± . (3-2) 

O o 

Bg 

Therefore this procedure generates a salinity profile from 
an observed profile consistent with the typical buoyancy 
profile, (3-1). Since this buoyancy profile is based on 
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the typical temperature and salinity profiles, it is repre- 
sentative of the conditions at that location during the 
specified month. 

The second method was to take the typical salinity pro- 
file and match it directly to the bathythermograph observa- 
tion. The advantage of the typical salinity profile over a 
climatological profile is that the typical profile is an 
actual observation and therefore has detailed structural 
features that the climatological profile is lacking. The 
mixed-layer depth of the BT was determined as above and the 
typical salinity profile was required to be mixed to the same 
or greater depth. The resulting profiles gave rise to a 
buoyancy profile which was then required to be stable: 

B(Z-AZ) >_ B(Z). If the profile was unstable at any level, 

Z, the salinity profile was adjusted by 

a [ T ( Z )-T (Z-AZ)] 

S T (Z) = S t (Z-AZ) + 2 ( 3- 3 ) 

which assures neutral stability. A new B(Z) was calculated 
and the procedure continued until the entire profile was 
stable . 

The salinity profile resulting from either of these 
methods permits positive temperature gradients that would 
appear unstable when considered alone. Including salinity 
therefore allows the model to be used over a much wider range 
of oceanic conditions without biasing potential energy as 
in Johnson (1977). 
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IV. DISCUSSION OF RESULTS 



Initially, the salinity flux at the surface of the ocean 
was neglected. This is potentially accurate only if there 
is a balance of evaporation and precipitation (E-P) over 
the time period of model integration. The validity of this 
assumption as well as the assumption of one-dimensionality 
can be checked by examining Figs. 1-3. In Fig. 1, the 
variability of surface salinity at OS PAPA between the 
first STP cast (1724 GMT 23 November 1974) and the last 
(1730 GMT 06 December 1974) is quite small and one can 
safely neglect the advection and surface flux of salinity. 
However, the variations of temperature and salinity in the 
pycnocline appear to be characteristic of advective affects 
associated with tidal period and internal wave motion. 

Also notice that the positive temperature gradients in the 
layer 80-150 m are excellent examples of the need for cor- 
responding salinity profiles when vertical mixing is being 
modeled. Variations at OS HOTEL between 1430 GMT 28 
November and 1442 GMT 29 December 1974 (Fig. 2) are on the 
order of one-half part per thousand and generally confined 
to the upper 50 m. Notice that there is more variability 
in the pycnocline here than at OS PAPA. At OS HOTEL in 
August (Fig. 3), salinity variations are quite extreme even 
within individual soundings. The variations cannot be 



25 



accounted for by realistic surface salinity fluxes and 
therefore are likely due to advective effects and/or small 
changes in ship position and the proximity to the Gulf 
Stream. 

A. HYDROCAST VS. BT INITIALIZATION 

The first model run shown using hydrocast data is at OS 
PAPA during the period from 1800 GMT 23 November to 1800 
GMT 06 December 1974. Wind stress and total heat flux. 

Figs. 4 and 5, respectively, for this period were extracted 
from the FNWC fields. The total heat flux is the sum of 
solar radiation, sensible and latent heat fluxes, and back 
radiation and is generally positive for this period, indi- 
cating an upward heat flux at the sea surface. Under these 
conditions, one would expect a deepening and cooling of the 
mixed layer. The initialization and verification profiles 
from STP data are shown in Fig. 6 and the predicted time 
evolution of surface salinity, surface temperature and mixed- 
layer depth are shown in Figs. 7, 8, and 9. The correspond- 
ing initial temperature profile (with positive gradients 
removed to ensure stability) derived from a BT and the pre- 
dicted changes in mixed-layer temperature and depth are 
depicted in Figs. 10, 11, and 12. 

Comparison of Figs. 8 and 11 shows a similar cooling 
trend and good agreement with the data. Notice that the 
initial SST from the STP cast was slightly warmer than the 
SST from the BT. Similarly, in Figs. 9 and 12, one can see 
that the STP initialization predicts a consistently 
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shallower MLD as compared to the BT initialization but each 
forecast agrees well if verified against similar type of 
data (see Figs. 6 and 10). This discrepancy in observed 
thermal structure is seen to be fairly consistent in that 
the sea-surface temperature averaged 0.26 °C warmer and the 
mixed-layer depth averaged 11.4 m shallower in the STP cast 
data than the BT data (Table II). Taking this into account, 
the model performance is very similar throughout the period 
as the STP initialization resulted in a cooling of the 
mixed layer by 0.8 °C and a deepening of 9.8 m whereas the 
BT initialization resulted in a cooling of 0.78 °C and a 
deepening of 7.1 m. Note that the slab depth (depth of the 
isothermal layer) in Figs. 9 and 12 retreats and deepens in 
a similar manner in response to weak stress and downward 
heat flux as described by Johnson (1977). This would be an 
expected feature since salinity, in the absence of sources 
and sinks, would have no effect on the shoaling of the iso- 
thermal layer. 

The major contribution of salinity in this case is the 
stabilizing effect on the positive temperature gradient 
below the thermocline. The salinity gradient that allows 
this to occur increases the density jump at the base of the 
mixed-layer depth and tends to inhibit deepening of the 
layer. This case further points out the importance of the 
initial profile in the subsequent model predictions, as was 
also shown by Johnson (1977). 

The next case considers a 31-day period (28 November to 
29 December) at OS HOTEL. The synoptic variability of stress 
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(Fig. 13) and the corresponding upward heat flux (Fig. 14) 
again indicate a deepening regime. The initial temperature 
profile from the Nansen cast in Fig. 15 and from the BT in 
Fig. 19 are essentially identical. The initial salinity pro- 
file (Fig. 15) shows a relative maximum at approximately 
100 m, which should act as a barrier to deepening. However, 
below this level, salinity decreases with depth and would 
contribute to vertical overturning if not stabilized by the 
temperature profile. A gradual increase of surface salinity 
is predicted (Fig. 16) as more saline water is entrained at 
the base of the layer. This variation appears to be smaller 
than the natural variability in the salinity measurements . 

Almost identical cooling rates of the surface of 3.7 °C 
during the period were forecast by hydrocast initialization 
(Fig. 17) and BT initialization (Fig. 20). Qualitatively 
this follows the trend indicated by the observations , but 
forecasts sea-surface temperatures that are on the order of 
1.5 to 2.0 °C too cold. Likewise, a very similar trend in 
mixed-layer deepening is forecast by the hydrocast initiali- 
zation (Fig. 18) and by the BT initialization (Fig. 21). 

Both forecasts tend to indicate too deep a MLD in the first 
14 days as compared to observations but then tend to corre- 
spond to the data in the later half of the period. As men- 
tioned previously, the slab depth responds identically in 
both cases. At the end of the period, the Nansen cast ini- 
tialization is approximately 7 m shallower than the BT 
initialization. This difference is due solely to the 
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salinity structure because the same forcing is used in each 
case and the initial temperature profiles were nearly 
identical . 

It may seem inconsistent that both mixed layers would 
have the same temperatures, since the deeper layer would 
have entrained more cold water from below. However, given 
the relatively weak temperature gradient and a difference 
in depth of only 7 m, the extra water entrained in the 
deeper layer has an insignificant cooling effect on a water 
column that is on the order of 100 m deep. 

B. TYPICAL PROFILE INITIALIZATION 

The lack of adequate temperature-salinity observations 
over large parts of the ocean was motivation for attempting 
to use the salinity profiles derived in Section III. Both 
procedures were tested in order to determine which, if either, 
was more useful. 

The typical buoyancy profile method for deriving a 
typical salinity profile. Equations (3-1) and (3-2) was 
tested first. This method, while providing initial profiles 
that were stable, gave salinity profiles that were not in 
general representative of those observed in the ocean. That 
is, unrealistic salinity variations were created below the 
mixed layer due to the constraint that buoyancy be conserved. 
This is shown schematically in Fig. 22. The typical buoyancy 
profile in Fig. 22a is mixed from to Z^ to correspond with 
the MLD of the BT. The salinity profile derived from the 
typical buoyancy profile is shown in Fig. 22b. In the mixed 
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layer, buoyancy is constant, thus the derived salinity will 
vary according to the temperature by Equation (3-2). At 
the base of the mixed layer, mixing of the buoyancy profile 
has created an artificially strong gradient between the 
level where mixing stopped, Z^, and the next layer below, Z^. 
Thus salinity increases to compensate for the buoyancy jump. 
Below level Z^, buoyancy has returned to a weak gradient 
but, this level is still in the thermocline of the observed 
BT. Thus salinity again decreases until the temperature 
gradient becomes weak. Below this, salinity again increases 
because the buoyancy gradient is greater than ag times the 
observed temperature gradient. Thus, the second approach of 
using a typical salinity profile and adjusting it to the 
bathythermograph was judged superior in general. However, 
notice this approach has the disadvantage of producing a 
larger than normal salinity jump at the base of the mixed 
layer if the typical salinity profile has to be mixed to a 
greater depth to correspond to the BT. Therefore, at deeper 
depths, (> 100 m), it might still be more appropriate to 
use the buoyancy approach, although this suggestion was not 
tested in this study. 

The model was initialized with a typical salinity profile 
and bathythermograph at OS PAPA and run for the same 13-day 
period as above. The results were cooling of the mixed 
layer by 0.75 °C and a deepening of the mixed layer by 5.17 m 
which agree relatively well with the results obtained by 
initialization with a BT alone. One might expect a typical 
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profile to work relatively well at this location because 
the overplots shown in Fig. 1 indicate very little variation 
of salinity. 

By the same token, one would not expect a typical pro- 
file to work as well at OS HOTEL given the variations of 
salinity in Figs. 2 and 3. This was tested for the same 31- 
day period as above. Temperature and corresponding derived 
salinity profiles are shown in Fig. 23. Notice that above 
100 m the typical salinity profile has a stronger gradient 
than the actual profile in Fig. 15. Both the initial (Fig. 
23) and the predicted (Fig. 24) mixed-layer salinities are 
consistently about 1 °/oo too fresh. However, this differ- 
ence in salinity has only a relatively small effect on the 
predicted mixed-layer temperature in Fig. 25, which is simi- 
lar to the trend in Figs. 17 and 20 with a mixed-layer cool- 
ing of 3.6 °C. Likewise the predicted mixed-layer depth in 
Fig. 26 displays a similar deepening trend as in Fig. 18 
and is only 5 m different by the end of the period compared 
to the model run using the Nansen cast initialization. Thus 
the model has made a good forecast even though the typical 
salinity profile is rather different then the observed pro- 
file in the upper levels. The noticeable difference between 
Fig. 26 and Figs. 18 and 21 is that the mixing, as indicated 
by the slab depth is strongly inhibited by the salinity jump 
at the base of the isothermal depth. Because of the 0.2 °C 
temperature difference between the slab depths and the mixed- 
layer depth definitions, this indicates a weaker temperature 
gradient below the well-mixed layer than is observed. 
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Of particular importance in Fig. 23 is the model pre- 
dicted development of a positive temperature gradient below 
the mixed layer. That is, the model mixed-layer tempera- 
ture is 13.78 °C and the temperature increases to 13.91 °C 
at 7.5 m below the isothermal layer. A positive gradient 
was observed in the verification BT. This type of thermal 
structure was also demonstrated by Miller (1976). 

C. SOURCES AND SINKS OF SALINITY 

In this section, the assumption of a long term balance 
between evaporation and precipitation was eased by including 
surface salinity flux in the model. The evaporation rate 
was derived from FNWC's Evaporative Heat Flux (EHF) field. 
The precipitation rate was obtained from 12-hour accumulated 
precipitation calculated by FNWC’s primitive equation model 
and was assumed to be evenly distributed over the 12 hours. 

Incorporating the salinity flux at the surface, the 
model was initialized with the temperature and salinity pro- 
files in Fig. 15. The resulting model profiles, after 31 
days, are shown in Fig. 27. The significant changes in the 
predicted salinity profile when compared to the salinity 
profile predicted neglecting (E-P) are the agreement of pre- 
dicted and observed mixed-layer salinities and the absence 
of a relative maximum at the base of the mixed layer (refer 
to Figs. 15 and 27). The evaporation rate minus the preci- 
pitation rate is shown in Fig. 28 and the corresponding sur- 
face salinity in Fig. 29. The surface salinities in Fig. 29 
show variations in response to (E-P) as opposed to the 
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surface salinities in Fig. 16 which show only an increase 
with time due to entrainment of more saline water at the 
base of the layer. Comparison of sea-surface temperatures 
[Fig. 17 without and Fig. 30 with ( E— P ) ] shows indistin- 
guishable forecasts. Thus sea-surface temperature changes 
due to inclusion of (E-P) are apparently not significant 
enough over a layer of this depth for the model to resolve. 
Finally, the slab depth responds differently under the in- 
fluence of (E-P) (Fig. 31) than when (E-P) was neglected 
(Fig. 18). The most significant changes caused by (E-P) 
are the enhanced retreats (decreases in slab depth) during 
days 3, 10, and 18 when precipitation was significant, and 
inhibited retreats during days 8, 13, 17, 22, 25, and 26 
when (E-P) apparently offset net downward surface heat flux. 

Of particular note is precipitation on day 18 which caused 
a retreat that corresponded to the observations. The subse- 
quent slow deepening of the slab depth represents the gra- 
dual mixing of the stable layer formed by precipitation. 

Typical profile initialization runs, both with and with- 
out (E-P), were made at other ocean station locations. Due 
to the sparsity of observations, no firm conclusions can be 
made as to the absolute verification of the forecasts. How- 
ever, a relative comparison can be made among results obtained 
from the different methods of initializing the model. These 
results are summarized in Tables III and IV. The length of 
the forecast period was determined by the continuous forcing 
function data and the longest period possible between initial 
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and verifying thermal structure data. The different model 
results shown in Table III for August indicate little 
variability in predicted SST or MLD due to salinity effects . 
This might be expected, in that August is a month of strong 
downward heat flux and temperature effects would dominate. 
Thus one might be cautious about including a typical salinity 
profile under these conditions since there is a possibility 
of degrading the model performance. Included in Table IV 
are the results obtained during a period of net upward heat 
flux. One would expect salinity effects to be more important 
in a deepening regime. The data in Table IV indicate, for 
the most part, that inclusion of salinity results in a warmer 
and shallower mixed layer than would be the case with den- 
sity as a function of temperature alone. Notice that, at 
OS INDIA (57N 20W) and OS MIKE (66N 02E), the MLD exceeded 
300 m (the bottom of the model). This might indicate that 
either the "typical" salinity profile used was not really 
representative of the area or that horizontal advection can- 
not be neglected here. However, in high latitudes in winter 
the temperature gradient is extremely weak. In this case, 
the slab depth would be more indicative of the depth to which 
the mixing processes were penetrating. At OS INDIA, the pre- 
dicted slab depths are 212.5 m (TYPPRO), 200.0 m (E-P), and 
272.0 m (BT). At OS MIKE, the respective values are 281.2 m, 
277.5 m, and 240.0 m. Although at OS MIKE the model appears 
to have run well on temperature only, the initial BT was cor- 
rected for a positive temperature gradient which would have 
biased the potential energy of the water column. 
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V. SUMMARY AND CONCLUSIONS 



It was the purpose of this study to test the effects of 
salinity in a mixed-layer model of the ocean thermal struc- 
ture. This is desirable in that the ocean thermal structure 
in many areas can be explained only if salinity is considered. 
However, this salinity data is not always available even at 
ocean weather stations. 

The problem was approached using two methods: (1) ini- 

tialization with hydrocast data, and (2) initialization with 
a BT and a salinity profile derived from historical data 
and adjusted to correspond to the BT. First, evaporation 
and precipitation were assumed to be in balance over the 
period of integration and the effects of salinity structure 
on the thermal structure were observed. Then, evaporation 
and precipitation were included to observe the effects of 
a surface salinity flux on the thermal structure. The in- 
corporation of observed salinity profiles into the model did 
not result in significant changes in forecasts of SST and 
MLD during the August period tests. Similarly, for the same 
conditions, inclusion of adjusted typical salinity profiles 
did not cause any distinct changes. However, the effects of 
salinity structure were noticeable in deepening regimes. 
Qualitatively, model performance was slightly improved by 
using hydrocast data. Use of adjusted typical salinity pro- 
files and BTs also resulted in slight improvement even when 
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the typical profile was somewhat different than the observed 
profile. The general effects of salinity structure were to 
produce a warmer, shallower mixed layer than results when 
density is considered a function of temperature only. 

Addition of the salinity flux due to imbalances of eva- 
poration and precipitation did not show pronounced effects 
on SST or MLD in the short term. However, a long term up- 
ward buoyancy flux, E>P, tends to offset the slower entrain- 
ment rate caused by the inclusion of salinity in the density 
profile. The short term effects of (E-P) were most evident 
in the response of the slab depth as compared to the response 
when only the surface heat flux was considered. The retreat 
of the isothermal layer is enhanced during precipitation and 
the subsequent deepening rate is slowed. Additionally, the 
tendency of this layer to shoal under downward heat flux may 
be inhibited by evaporation. 

The addition of salinity effects to the model is certain- 
ly desirable. However, salinity data are sparse and the 
typical salinity profiles must be applied with an understand- 
ing of the meteorological and oceanographical conditions of 
the area. Thus one would consider salinity data to be far 
more important in mid and high latitudes during a deepening 
regime than in low latitudes during strong heating. One must 
also consider the variations of salinity and the processes 
causing these variations. In subpolar regions, salinity 
varies markedly with freezing and thawing. In monsoonal 
areas, reversal of the wind and thus wind-induced currents, 
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makes advective effects important. Further, an attempt to 
apply the model ocean-wide would require salinity data in 
many areas where it is not now available. 
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APPENDIX A 



HYDROCLIMATOLOGICAL DATA 
RETRIEVAL PROGRAM 

The Hydroclimatological Data Retrieval (HYDAT) program 
was designed by Ocean Data Systems, Inc. (1977) to access, 
process, display, and store hydroclimatological data bases. 
Among the outputs available from HYDAT is an option that 
produces a typical temperature profile. This typical profile 
is obtained by specifying the month (or months) and the loca- 
tion of interest. The location specified is taken to be the 
midpoint of a circular or square search area, with the 
original search distance from the midpoint determined by the 
user. If at least six reports are not found in this area, 
the area is expanded by fixed increments until six or more 
reports are found or until the maximum search distance, 
either default or specified, is reached. 

If the search is successful, each report found is sepa- 
rated into three parameters: (1) the sea-surface temperature, 

(2) the deep temperature, a temperature at an arbitrary depth 
('•2 00 m) , and (3) the heat content as computed from 

(A-l ) 

i=l 

where D is the deepest level, N is the number of layers to 
depth, D,AZ is the thickness of each layer, and T is the 
mean temperature in the layer. Any sounding not extending 
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to D is not considered in the selection of the typical pro- 
file. The values of each • parameter are then screened to 
remove isolated extreme values and the frequency distribu- 
tion is checked for the number of modes (£3) present. All 
three parameters must have the same number of modes. HYDAT 
then computes and prints the median and first and third 
quartile values for each mode of each parameter. The 
typical profile in each mode is determined by comparing the 
parameter values of each profile with the median values and 
determining a deviation, D . , 



D. 

1 





i=l,2,3 



(A- 2 ) 



where is a weighting factor assigned to the individual 
parameter, is the value of the parameter, and is the 
median value of the parameter within any one mode. The 

typical profile is then the report with the lowest sum of 

/ 

deviations : 



Sum = 2 D . 

l 



l 



( A- 3 ) 



TABLE I. Data availability at ocean station locations. 

Numbers under TYPICAL PROFILE indicate the 
total number of profiles (both BT and hydro- 
cast) upon which the typical profile is based. 
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TABLE II. Comparison of STP casts with BTs at OS PAPA, 
November-December 1974. Observations are 
within one hour and ten nautical miles of 
each other. 





SS rdrp) 


s 5 rear) 


' AT 


/m-DCstP) 


M>-0(feT) 


AD 


2.3 Me v 


8.11 °C 


8 .D e C. 


6 . il # C. 


♦ 

6>5". 3 m 


7k. 2 . no 


-10.7 m 


^4 Mcv/ 


7.40 


7.k 


0.3 


20 . 4 * 


73. 8 


• 2 , c l 


2 k Dc y 


7. 87 


7.t 


0 . 2.1 


7 Z.. 1 


6fe.7 


— Ik, (o 




7. C H 


Lto 


©31 


11.0 


8 3.4 


-12.4 


01 bee 


7.84 


7.(o 


0.24 


kb, 1 


61.3 


-i(p.2_ 


62.0to 


1.^3 


7 7 


o. z.3 


12.3 


11.0 


o.3 


o3De<_ 


7.<ol 


7.4 


0 . 2.1 


13" .4 


m3*" 


- 12,4 


C>4 Oee 


7.8| 


IS 


o.3J 


kb. 4 


&ST.<j 


-io.S' 


Ot Dec 


7 -k 6 


1M 


0.28 


Ife.k 


8 t.&' 


-)0,’L 



fAoxv\ O.Zk C~ HMff) 



^Tndicde*, pciSi'/u/c cifGcii(^o4 c\l6ove M*-D, 



41 



TABLE III. 



Comparison of forecast SST and MLD among 
different model initialization runs and 
with observations for August 1974. HC is 
hydrocast data, TYPPRO is BT and typical 
salinity data, (E-P) indicates inclusion 
of surface salinity flux, and BT is tem- 
perature data. Verifying data are from BTs . 
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TABLE IV. Comparison of forecast SST and MLD among 
different model initialization runs and 
with observations for November-December 1974 . 
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Figure 1. Overplot of hydrocast data for OS PAPA for 
the period 21 November-31 December 1974. 
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Figure 2 . Overplot of hydrocast data for OS HOTEL for 
the period 21 November-31 December 1974. 
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Figure 3. Overplot of hydrocast data for OS HOTEL for 
the period 8-28 August 1974. 
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Figure 4. Wind stress at OS PAPA - 18 GMT 23 November 
to 18 GMT 06 December 1974. 
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Figure 5. Total heat flux at OS PAPA - 18 GMT 23 November 
to 18 GMT 06 December 1974. SLB is the sum of 
sensible and latent heat and back radiation. 

THF is SLB plus solar radiation. 
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Figure 6. Depth vs temperature and salinity at OS PAPA - 
18 GMT 23 November and 18 GMT 06 December 1974, 
hydrocast data, neglecting (E-P). INITIAL in- 
dicates the initialization profiles, MODEL in- 
dicates the models predicted profiles , and OBS 
indicates the observed profiles at the end of 
the period. 



49 



33.30 




2 




(2 


r s 


(2 


12 


(2 


2 


(2 


CD 




cs» 


S3 


CO 


CD 




CM 


t2 


* 


• 


• 


• 


% 


9 


% 


• 




on 


m 


on 


CO 


CM 


CM 


CM 


CM 


CM 


CO 


CO 


CO 


00 


CO 


00 


CO 


CO 


CO 



-st - 

f>- 



O 



M 

CO 

r—H 

CO 

CM 



<3 

cn 



CO 

'jt — 

CE 

O 



(idd) AlINIldS 



50 



Figure 7. Surface salinity at OS PAPA - 18 GMT 23 November to 18 GMT 06 

December 1974, neglecting (E-P). MODEL is the model prediction 
of salinity and OBS indicate observed values. 
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Figure 8. Mixed-layer temperature at OS PAPA - 18 GMT 23 November to 18 
GMT 06 December 1974 from hydrocast initialization. MODEL in 
dicates predicted temperature, OBS indicate observed values, 
and PG indicates an observed positive temperature gradient. 
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Figure 9. Mixed-layer depth at OS PAPA - 18 GMT 23 November to 18 GMT 
06 December 1974 from hydrocast initialization. SLAB is the 
model depth of the isothermal layer, MODEL is the model MLD , 
and OBS indicates observed values . 
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Figure 10. 



Depth vs temperature at OS PAPA - 18 GMT 23 
November and 18 GMT 06 December 1974, BT data. 
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Figure 11. Mixed-layer temperature at OS PAPA - 18 GMT 23 November to 
18 GMT 06 December 1974 from BT initialization. 
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Figure 12. Mixed-layer depth at OS PAPA - 18 GMT 23 November to 18 GMT 
06 December 1974 from BT initialization. 
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Figure 13. Wind stress at OS HOTEL - 14 GMT 28 November to 14 GMT 
29 December 1974. 
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Figure 14. Total heat flux at OS HOTEL - 14 GMT 28 November to 
14 GMT 29 December 1974. 
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Figure 15. Depth vs temperature and salinity at OS HOTEL 
- 14 GMT 28 November and 14 GMT 29 December 
1974, hydrocast data, neglecting (E-P). 
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Figure 16. Surface salinity at OS HOTEL - 14 GMT 28 November to 
14 GMT 29 December 1974, neglecting (E-P). 
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Figure 17. Mixed-layer temperature at OS HOTEL - 14 GMT 28 
November to 14 GMT 29 December 1974 from hydro- 
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Figure 18. Mixed-layer depth at OS HOTEL - 14 GMT 28 November to 
14 GMT 29 December 1974 from hydrocast initialization. 
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Figure 19. Depth vs temperature at OS HOTEL - 14 GMT 28 
November and 14 GMT 29 December, BT data. 
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Figure 20. Mixed-layer temperature at OS HOTEL - 14 GMT 28 November 
to 14 GMT 29 December 1974 from BT initialization. 
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Figure 21. Mixed-layer depth at OS HOTEL - 14 GMT 28 November to 
14 GMT 29 December 1974 from BT initialization. 
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Figure 22. Schematic of a derived salinity profile 
assuming a typical buoyancy. 
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Figure 23. Depth vs temperature and salinity at OS HOTEL 
- 14 GMT 28 November and 14 GMT 29 December 
1974, BT and typical salinity data, neglecting 
(E-P). 
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Figure 24. Surface salinity at OS HOTEL - 14 GMT 28 November to 
14 GMT 29 December 1974, neglecting (E-P) and using 
a typical salinity profile. 
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Figure 25. Mixed-layer temperature at OS HOTEL - 14 GMT 28 
to 14 GMT 29 December 1974 from typical profile 
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Figure 26. Mixed-layer depth at OS HOTEL - 14 GMT 28 November to 
14 GMT 29 December 1974 from typical profile initiali- 
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Figure 27. Depth vs temperature ' and salinity at OS HOTEL 
- 14 GMT 28 November and 14 GMT 29 December 
1974, hydrocast data. 
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Figure 28.* (E-P) at OS HOTEL - 14 GMT 28 November to 14 GMT 

29 December 1974. 
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Figure 29. Surface salinity at OS HOTEL - 14 GMT 28 November to 
14 GMT 29 December 1974. 
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Figure 30. Mixed- layer temperature at OS HOTEL - 14 GMT 28 November 
to 14 GMT 29 December 1974. 
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Figure 31. Mixed-layer depth at OS HOTEL - 14 GMT 28 November to 
14 GMT 29 December 1974. 
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